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The c-axis optical response of the underdoped cuprates is qualitatively different from its in-plane 
counterpart. The features of the pseudogap show themselves more prominently in the c-axis than 
in-plane. We compute both the c-axis and in-plane optical conductivity using the Yang-Rice-Zhang 
model of the underdoped cuprates. This model combined with coherent interlayer tunnelling is 
enough to explain the qualitative differences between the in-plane and c-ax;is data. We show how 
pseudogap features manifest themselves in the infrared and microwave conductivity within this 
model. 
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I. INTRODUCTION 

The nature of the pseudogap phase in the underdoped 
cuprates is believed to be central to the understanding 
of high- Tc superconductivity.— Many ideas have emerged 
to understand the origin of the pseudogap. Examples 
of theories include the idea of preformed cooper pairs^ 
or a competing order parameter, such as d— density wave 
order. An alternative picture has its roots in Anderson's 
resonating valence bond (RVB) order. ^'^ Within the RVB 
framework the pseudogap can emerge naturally as one 
dopes a Mott insulating state with holesi^ The model 
of the underdoped cuprates by Yang, Rice, and Zhang 
(YRZ) is based on these ideas.'' 

Since the YRZ model was put forward, it has proved 
successful in describing many features of the underdoped 
cuprates that cannot be understood from conventional 
BCS theory. The essential new feature is the pres- 
ence of an additional energy scale, namely, the pseu- 
dogap. In the YRZ model the pseudogap is responsi- 
ble for reconstructing the antinodal portion of the Fermi 
surface into closed Luttinger pockets. With this modi- 
fication, the YRZ model has been able to qualitatively 
capture the physics of Raman Spectraj^i^ ARPESfiSiiii 
specific heatr^ penetration depthfi2. as well as tunneling 
spectroscopy^^ More recently, we have applied the YRZ 
model to the c-axis transport properties where we showed 
that the YRZ model is able to explain the insulting-like 
c-axis behaviour while remaining metallic in-plancjA'^ It 
has also been shown to account for the c-axis violation of 
the Ferrell-Glover-Tinkham sum ruleJ^ It is remarkable 
that such a simple modification is capable of capturing 



the physics of such a diverse range of topics. The recent 
ARPES observation of closed pocketsii adds further sup- 
port for the YRZ model. 

The optical response of the c-axis is known to be dra- 
matically different from the in-plane response, both in 
for infraredi^ as well as for microwave^i^ frequencies. In 
this paper we examine the differences between in-plane 
and c-axis optical response. We compute the optical con- 
ductivity using the YRZ formalism. We show that the 
YRZ model is able to capture the qualitative behavior of 
the AC optical conductivity both in-plane, and along the 
c-axis. We use the conductivity to extract information 
about the behaviour of the superfluid density, as well as 
the distribution of optical spectral weight. Lastly, we 
use the low frequency portion of our data to extract the 
microwave conductivity and find good agreement with 
experimental findings. In section |ll] we introduce the 
formalism required to compute the optical conductivity 
within the YRZ model. We present our numerical re- 
sults for the infrared optical conductivity in section IIIIl 
and discuss optical sums and the microwave conductivity 
in section IIVI We summarize and conclude in section |Vl 



II. OPTICAL CONDUCTIVITY IN THE YRZ 
MODEL OF THE UNDERDOPED CUPRATES 



The real part of the c-axis optical conductivity in 
the bubble approximation can be expressed in terms of 
the spectral density, A{\<i,uj), and the Gorkov anomalous 
spectral density, i3(k, w), through the Kubo formula: 



Re[cj,{uj,T)] 



duj' 
"2^ 



[/(w' + oj)~ f{uj')] [A{k, uj')A{k, uj' + uj) + B{k, uj')B{k, uj' + uj)] . (1) 
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Here e is the electron charge, /(w) is the Fermi distri- 
bution function, t_L(k) is an interlayer hopping matrix 
element, and d is the interlayer distance. We will often 
make comparisons with the in plane conductivity where 
t'^d'^ should be replaced by vl , the electron velocity. The 
YRZ model provides the coherent part of the Greens 
function from which we can extract the required spec- 
tral densities. For a doping, x, the Greens function is 
given by 



(2) 



In the above gt {x) is a Gutzwiller renormalization factor 
and is given by gt{x) = 2x/(\ + x). The two energy 
branches and weights are given by 
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= o(ek-ek)±i?k, 
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-2t{x){coskx + cos ky), 
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~ 4i'(x) cos kx cos ky 

— 2t" (x) (cos 2kx + cos 2ky 



Hp. 
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Hp is a chemical potential that is determined from 
the Luttinger sum rule. The hopping parameters are 
t{x) = gt{x)to + {3/8)gs{x)Jx, t'{x) = gt{x)t'o, and 
t"{x) = gt{x)tQ, where gs{x) = 4/(1 + x)"^ is another 
Gutzwiller renormalization factor, J = to/3, x = 0.338, 



-O.Sio and t'^ = 0.2io. 



In Eq. ([S|) Apg is the pseudogap energy scale which 
is taken to have d-wave symmetry, along with the super- 
conducting gap. That is. 



A.,. - ""^ ^ (cosfc^ 



2 



(cos kx 



cos ky), 
cos ky). 



(9) 
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The doping dependent magnitudes mimic a simplified 
version of the cuprate phase diagram (see left panel of 
Fig. [J): 

A°^(a;) =0.14(1 -82. 6(x- 0.2)2) ^^^y^ 
A°g(a:)= 0.6(1 -x/0.2). (12) 




FIG. 1: (Color online) The phase diagram and Fermi surface 
for the Yang-Rice-Zhang model. The quantum critical point 
at which the pseudogap emerges is set to 2; = 0.2 in this work, 
and corresponds to optimally doped superconductivity. As 
the pseudogap grows with reduced doping the Fermi surface 
reconstructs from the large band structure Fermi surface at 
optimal doping into closed Luttinger pockets. 



This functional form places both optimal doping (the 
maximum of the superconducting dome, shown in dashed 
purple in Fig. [Ij and the vanishing of the pseudogap en- 
ergy scale (shown by the solid blue line in Fig. [1]) at 
X = 0.2, in accordance with the original YRZ paper—. 
Unless otherwise specified, we use all parameters in the 
band structure for the YRZ model as they appear in the 
original publication.^ In principle, one could alter these 
parameters to obtain fits to experimental data, but in this 
work we only wish to show that YRZ captures the essen- 
tial physics. For the magnitude of the superconducting 
gap, we use the ratio 2A°^(a;,T = 0)/(fci3T'c) = 6. 

From the Greens function we obtain the spectral den- 
sities 



A(k, uj) = 2i:g^(^x) J] {in'^fSiuj ~ E'^) (13) 

+(z;")25(c. + £;,")], 



a=± 



and. 



B{\^,uj) = 2^gt{x) W^iu^v^ [6{^ ~ Ef) - 6{u; + Ef)] 

(14) 



a=± 



The Bogoliubov quasiparticle energies, Ef, and ampli- 
tudes, u" and are given by 

(15) 
(16) 
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FIG. 2: (Color online) (Top panel): The optical conductivity 
at doping x — 0.13 and T = 0.25Tc. Both the ab-plane and 
c-axis contain a coherent Drude response (red) and an inter- 
band contribution from the presence of the pseudogap energy 
(orange). (Bottom panel): The optical conductivity at opti- 
mal doping (x = 0.20) and T = 0.25rc. At optimal doping the 
pseudogap vanishes and there is only a Drude contribution to 
the conductivity. 

The right hand frame of Fig. [1] shows the normal 
state (A°^(a;) = 0) Fermi surface reconstruction brought 
about by the emergence of the pseudogap. Only the up- 
per right quadrant of the Brillouin Zone (BZ) is shown. 
The dashed black line indicates the antiferromagnetic BZ 
(AFBZ) boundary. The solid green curve in the upper 



frame is the large Fermi surface of Fermi liquid theory for 
a doping of x — 0.2. The Fermi liquid contour crosses the 
AFBZ but is unaffected by it and is characteristic of a 
good metal. As the doping is reduced towards half filling 
the Mott insulating state is approached and the Fermi 
surface contours change radically. The lower frame is 
for a doping of x = 0.13 where the Fermi surface has 
reconstructed into a Luttinger hole pocket. The back- 
side of the Luttinger pocket at a; = 0.13 closely follows 
the AFBZ boundary and has very small weight {Wk in 
Eqn. (U])). On the other hand, the front side is heavily 
weighted and is very close to the underlying Fermi liquid 
surface (Apg = 0) in the nodal direction. As x is reduced 
further the Luttinger pocket continues to shrink and only 
a small number of well defined quasiparticles remain in 
the nodal direction. This fact is very important for much 
of the physics that we will describe in this paper. 

These small hole-like pockets are a prediction of the 
YRZ model and are in excellent agreement with recent 
photoemission datafii However, it is believed that some 
of the transport properties in the pseudogap phase are 
electron-like. Experimentally, it is found that both the 
Hall and Seebeck coefficients are negative in the pseudo- 
gap phase. This apparent contradiction is nicely over- 
come in a recent proposal by Harrison and Sebastian.— 
They show how a nodal hole-like Fermi arc can be 
reconstructed into electron-like pockets. Their recon- 
struction mechanism is due to bilayer charge ordering. 
The wavevectors associated with the charge modulation 
are responsible for the reconstruction from the hole-like 
Fermi surface into the electron-like one. Importantly, the 
reconstruction happens with the nodal piece of Fermi sur- 
face in their model. The YRZ model naturally makes hole 
pockets with long lived quasiparticles along the nodal 
direction. Additionally, YRZ predicts that these hole 
pockets are heavily weighted along the front side of the 
pocket, and so this reconstruction mechanism should ap- 
ply in exactly the same way as they describe. This ex- 
tra reconstruction offers a nice explanation for why these 
transport coefficients are observed as negative in this part 
of the phase diagram. 

Using the above spectral functions, the conductivity 
can be written as the sum of two terms, cr = gu + crjB- 
The first term, an, is peaked around uj — and is a 
Drude-like response, while the second term, ajB arises 
from interband transitions between the different energy 
branches. In our calculations we take into account the 
effect of impurities by replacing the Dirac delta functions 
by Lorentzians of half-width F. In the clean limit we find 



-2 df{E7) ' 
dE7 



(18) 
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(19) 



When using the clean hmit formulas, one must instead 
replace the Dirac delta functions by Lorentzians of half- 
width 2T. 

For the interlayer tunneling matrix element we follow 
Anderson^" and choose t±{'k) = t± {cos{kx) — cos(fcj^))^. 
This choice reflects the geometric arrangement of the 
atoms between adjacent Cu02 planes. For the ab-plane 
the velocity is simply Vk^ = di^ldkx. The last free pa- 
rameter we have is the scattering rate which broadens 
the Dirac delta functions. We used 

fo.OOl + O.l (^)' T<T„ ^20) 
[0.051 + 0.05^ T>Tc. 

A linear in T quasiparticle scattering rate has been asso- 
ciated with the inelastic scattering in the high oxides. 
It is taken to be a characteristic of their normal state. 
The marginal Fermi liquid phenomenology^! is based on 
this observation as well as the idea that the the dominant 
scattering processes involve spin and charge excitations 
of the electronic system itself. On entering the supercon- 
ducting state, the emergence of the superconducting gap 
reduces both the charge and spin susceptibility and leads 
to a reduction in scattering. This reduction in scattering 
is a hallmark of an electronic mechanism for the inelastic 
scattering and is often referred to as the collapse of the 
inelastic scattering ratei^"— This collapse of the scat- 
tering rate is the accepted explanation of the large peaks 
observed in the microwave^l and thermal conductivity^^ 
of the cuprates well below T^. While the normal fluid 
density (which is resistive) is reduced with temperature, 
the inelastic scattering lifetime increases. The increase 
ceases when the residual scattering becomes dominant 
and it is the further reduction in normal fluid density 
that drives the conductivity to zero as observed in the 
experiment 1^ 



III. NUMERICAL RESULTS 

After all of these choices we can evaluate the conductiv- 
ity. Fig. [2]shows the result of a calculation at T = 0.25rc 
for optimal doping {x — 0.20) as well as underdoped 
(x — 0.13). We chose x = 0.2 to highlight the physics in 
the absence of a pesudogap. By contrast, the x = 0.13 
case has sizable pseudogap and superconducting energy 
scales. This choice allows us to see the effect of both en- 
ergy scales in our data most easily. In the x = 0.20 case 



(bottom two frames of Fig. (2) , there is only a coherent 
Drude-like response from the large Fermi surface in both 
the ab-plane and the c-axis. The c-axis response is re- 
duced as compared with the ab response by the out of 
plane matrix element (k) which gives less weight to the 
part of the Fermi surface in the nodal direction. In the 
underdoped case (top two frames of Fig. [5]), there is still 
a Drude response from the remaining Fermi surface. It is 
suppressed in both the ab-plane and c-axis due to both 
the shrinking size of Fermi surface and the Gutzwiller 
renormalization factors. There is also a piece due to in- 
terband transitions at higher energies. This piece is the 
signature of the pseudogap energy scale. It is the domi- 
nant feature in the c-axis data. 

In more realistic treatments of the inelastic scattering 
we expect this interband feature to broaden beyond what 
is shown in Fig. [51 In the cuprates the inelastic scatter- 
ing rate is known to have a strong energy dependence— 
in addition to temperature variation. In Eqn. (j20p we 
have included only temperature dependence through an 
impurity model, i.e., there is no dependence on uj. This 
model for the inelastic scattering rate is perfectly ade- 
quate as far as DC properties are concerned. It is over- 
simplified when the photon energy falls in the infrared. 
In this energy range, T{T,uj) can be much larger than 
its w = value, modeled in Eqn. ((20)) . Consequently, 
the interband optical transition peak in Re[cr(t,aj)] will 
be broadened. 

In Fig. [3] the optical conductivity is shown as a func- 
tion of frequency for four temperatures: as labeled in 
the figure, T = 4Tc dashed blue, T = 2Tc solid orange, 
T = Tc in solid green with boxes and T — 0.25Tc in 
solid black. The top two frames are for doping x — 0.13 
and the middle two at optimal doping (x = 0.2). For 
Lo < 0.2^0 all curves are Drude like and order in temper- 
ature in the usual way, except the c-axis curves in the 
underdoped case. In this case, we see a different trend, 
the curves order in the opposite sense as temperature 
is reduced. This behaviour is emphasized in the lowest 
frame where we show our results (for x = 0.05) for the DC 
resistivity as a function of temperature. We see that the 
c-axis displays insulating like behavior {pc increases as T 
decreases), while the ab-plane response remains metallic 
{pab decreases as T decreases). 

The top two frames of Fig. [3] arc to be compared with 
the experimental work of Homes et al. on the c-axis 
conductivity of YBCO.^^ As in the experiments, the sig- 
natures of the pseudogap are more pronounced at low 
energies. This is characterized by a fiat region at low fre- 
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FIG. 3: (Color online) (Top panel): The optical conductivity 
at doping x = 0.13 for temperatures T — 0.25Tc, Tc, 2Tc, and 
4:Tc. The Drude peak becomes increasingly broadened with 
temperature. At low temperatures a feature of the pseudogap 
emerges beyond u — O.Sto- This feature is a natural expla- 
nation for the broad peak observed at 400cm~^ in Homes 
et al.r"^ (Middle panel): The optical conductivity at dop- 
ing X = 0.20 for temperatures T = Q.2hTc,Tc,2Tc, and ATc. 
Both in-plane and c-axis only have Drude contributions which 
become increasingly broadened with temperature. (Bottom 
panel): The resistivity in the underdoped cuprates from the 
YRZ model as a function of temperature at a: = 0.05. The 
underdoped phase is metallic in-plane and resistive along the 
c-axis. Taken with permission from Ashby and Carbotte.— 
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FIG. 4: (Color online) Normal state optical sum minus su- 
perconducting state optical sum as a function of temperature 
for doping x — 0.13. The behaviour of the superfluid density 
in the c direction is qualitatively different from the ab-plane. 



quency which becomes increasingly suppressed at lower 
temperatures. The reduction we calculate is not as dra- 
matic as that observed in the experimental data. If one 
wished to obtain better fits to the experimental results it 
could be achieved by adjusting the scattering rate or the 
magnitude of the pesudogap energy scale. Note that the 
low frequency c-axis conductivity behaves very differently 
from the in-plane conductivity as a function of temper- 
ature. Perhaps the most striking feature in Fig. |3] is 
the broad interband feature from lj ^ O.Mq to a; ^ O.Stg 
associated with the pseudogap. This feature is a natu- 
ral explanation for the broad peak observed in Homes et 
ali^ at 400cm^^. Taking to — 125meV, this feature falls 
in the range 300 — 500cm^^, in agreement with the ex- 
periments. Previous explanations of this feature invoke 
the existence of an interlayer plasmon collective modei^ 
Within YRZ this peak is produced only from pseudogap 
physics with no need for this collective mode. 



IV. OPTICAL SUM AND MICROWAVE 
CONDUCTIVITY 

A. Optical sums 

Sum rules and partial optical sums also provide use- 
ful information. In a conventional superconductor, the 
suppression of Re[cr] for T < Tc is connected to the ap- 
pearance of a superconducting condensate. The 'missing' 
optical spectral weight appears in a delta function of the 
appropriate weight at a; = 0. This fact is usually pre- 
sented as the Ferrel-Glover-Tinkham sum rule. In terms 
of the superfluid stiffness, p^: 



N 



N 



(21) 
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where 



iV^v - / dcjRe[(7„(w,T)] 
is the normal state optical sum, and 



0+ 



(iwRe[(Ts(cj, T)] 



(22) 



(23) 



is the superconducting optical sum. While this sum rule 
holds in-plane, it is known to be violated in the c-axis<^ 
The YRZ model with the same interlayer tunneling ma- 
trix element that we use displays a violation of this sum 
rule This is understood most simply in a limit of YRZ 
that reduces to a Fermi arc model. In this model the 
Fermi surface is confined about the nodal direction; the 
remainder of the large Fermi surface from Fermi liquid 
theory is gapped out by the pseudogap. Only the elec- 
trons on the arc contribute the usual amount to the op- 
tical sum rule, while those on the gapped out portion 
contribute less. In any case, — Ns is well defined 
and is shown in Fig. The solid blue line is the ab- 
plane result, while the orange dashed line applies to the 
c-direction. As a function of the reduced temperature, 
t = T/Tf. those curves behave in much the same way 
is is found for the superfluid density itself The su- 
perfluid density follows from the imaginary part of the 
conductivity as 



1 



47r 



A2(T) 



— lim — a;Im[(T(T, cj)]. 



>o c 



(24) 



with c the velocity of light. The temperature behaviours 
are in agreement with the superfluid density inferred from 
microwave experiments The ab-plane superfluid den- 
sity decreases linearly with temperature while the c-axis 
is flat at low temperatures. 

A related quantity, the partial optical sum 



0+ 



dw' (Re[(T„(c^',r)] - Re[(T,(u;',r)]) , (25) 



provides information about the distribution of spectral 
weight that goes into the superconducting condensate. 
The partial optical sums for the YRZ model, normalized 
to their value at w = 1 are shown in the top panel of 
Fig. [5] The solid blue curve is for the ab-plane and 
the dashed orange curve is for the c-axis. For a regular 
superconductor this curve would sharply rise towards 1, 
with the scattering rate, F, setting the energy scale of 
the rise. In the YRZ model the presence of the pseudo- 
gap causes a redistribution of spectral weight from the 
Drude to the high frequency region. This redistribution 
shows itself in the suppressed flat region at low frequen- 
cies accompanied by a pile up of spectral weight above 
the pseudogap energy scale. This redistribution is much 
more pronounced for the c-axis. Insight into these shifts 
in spectral weight can be gleaned from the lower frame 
where we show the real part of the conductivity as a func- 
tion of frequency for x = 0.13 at T = 0.25rc- The dashed 
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FIG. 5: (Color online) (Top Panel): Partial optical sum at 
T = 0.25rc for doping x — 0.13 for the ab-plane and the 
c-axis. The gap scales are indicated on the Fig.. The pseudo- 
gap causes an extra suppression and redistributes the charge 
carriers available for condensation to high frequencies. (Bot- 
tom Panel): c-axis conductivities in the superconducting and 
normal states for s = 0.13 at T = 0.25Tc. The redistribution 
of spectral weight due from having superconductivity leads to 
a shift in the interband peak. 



orange curve is the normal state (Ap^ = 0) and the black 
curve is the corresponding superconducting case. We see 
that the opening of the superconducting gap shifts the 
interband transitions to higher energies. This shift ac- 
counts for the large peak seen in the orange curve in the 
top frame of Fig. [5]for > QAto. 



B. Microwave conductivity 

The microwave data is qualitatively different for the 
in-plane and c-axiSi^S The in-plane data contains a peak. 
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FIG. 6: (Color online) Microwave conductivity for a constant 
r = 0.005 as a function of temperature. This is proportional 
to the number particles able to absorb low frequency radia- 
tion. The behavior in plane is linear in T, while the c-axis 
behaves roughly like . 



while the c-axis does not. The peak is attributed to 
the interplay between the quasiparticle lifetime and the 
amount of fluid which can absorb radiation. As the sys- 
tem is cooled below Tc the quasiparticle scattering life- 
time changes. As the temperature decreases the nor- 
mal component of the fluid vanishes linearly (Fig. [S]). 
If the lifetime increases faster than the normal part of 
the fluid, this will lead to an increasing conductivity. At 
low enough temperature, the lifetime will saturate to the 
value set by the residual scattering rate, and the conduc- 
tivity will decrease as the remaining normal fluid con- 
denses. This behavior generally produces a peak in the 
microwave data. It is thus, very surprising that the c-axis 
lacks a peak. This was originally interpreted as evidence 
for incoherent c-axis transporti^. In this work we inter- 
pret the difference as a signature of how the superfluid 
density is changing. 

The microwave conductivity is obtained from the low- 
frequency part of Re[CT(a;)]. We obtain the microwave 
conductivity by taking lim^j-^o Re[tT(cij)]. In Fig. |6] we 
show the microwave conductivity for a constant scatter- 
ing rate F = 0.005. The solid black curve is Re[cr(T, a; = 
0)] for the ab-plane while orange is for the c-axis both as 
a function of reduced temperature, T/Tc- Their temper- 
ature dependence is striking different as we would have 
expected based on the results presented in Fig. |4] for 
Niq — Ns vs T/Tc- For a first understanding of the ex- 
perimental results for Re[i7(T, w = 0)], it is helpful to 
take guidance from the two fluid model. Under the as- 
sumption that only the normal fluid component, n, is 
involved in the absorption, the expression for the Drude 



FIG. 7: (Color online) The in-plane and c-axis microwave 
conductivity for x — 0.20 as a function of temperature for 
various different forms of scattering rate. 
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FIG. 8: (Color online) The microwave conductivity in the 
underdoped phase {x — 0.13) as a function of temperature. 
We also show the calculation with Apg = 0. The pseudogap 
has almost no effect on the in-plane conductivity, but leads 
to a dramatic suppression for the c-axis. 



conductivity is 



mV 



(26) 



Here e is the electron charge and m is the mass. We see 
in Fig. |6]that n is nearly linear in plane and oc out of 
plane. 

In Fig. [7] we show the microwave conductivity at 
optimal doping for different scattering rates. We used 
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r = 0.001 + 0.1(T/Tc)" for a = 2,3,4. Dots correspond 
to law, dashed to and solid to T^, with black for ab 
and orange for c-axis. The in-plane data always contains 
a peak, as we expect. If the scattering rate becomes too 
strongly temperature dependent, the c-axis conductivity 
is no longer convex like in the experiments. Using the 
Drude form for the conductivity with T = A-\- BT"^ , and 
n oc we can show that for the conductivity to be con- 
vex that a < 1 - 2^ + 2v/2/3(/3 - 1). Using /3 w 5 from 
our constant scattering rate calculation, we see this sets 
an upper bound a = 3.65. 

In Fig. |5] we show the microwave conductivity for 
X = 0.13. To isolate the effect of the pseudogap, we 
redid the calculation with Ap^ = 0. The solid curves are 
with the pseudogap and the dashed without. The low 
temperature behavior is unaffected by pseudogap forma- 
tion. It does cause a dramatic suppression of the c-axis 
conductivity at high temperatures, but this effect is much 
more modest in-plane (Fig. There is a great similar- 
ity between the change in temperature behavior of the 
c-axis microwave conductivity and the specific heatj^ In 
both cases the low temperature part of the curve is un- 
affected since this region depends only on the thermal 
excitations in the nodal direction. This part of the elec- 
tronic structure is not appreciably changed by the pseu- 
dogap. As the temperature is increased towards the 
specific heat is strongly suppressed below its Apg = 
value, much like the c-axis conductivity. This is not sur- 
prising as both quantities are closely tied to the electronic 
density of states. The in-plane microwave conductivity 
is not, and is seen to behave much differently from its 
c-axis counterpart. Unfortunately c-axis measurements 
are technically challenging and experimental data only 
exists for optimal doping. It would be very interesting 
to look for the effect of the pseudogap in an underdoped 
sample. 

V. DISCUSSION AND CONCLUSIONS 

We have investigated the c-axis optical conductivity in 
the underdoped cuprates using the YRZ model. We fo- 
cused on properties in the superconducting phase of the 
underdoped cuprates at x — 0.2 (optimal doping) and 
X = 0.13 (underdoped) to highlight the essential features 



of the model. For the c-axis calculations we used a co- 
herent tunneling matrix element to describe interlayer 
hoping. Our choice of matrix element is one related to 
the geometric alignment of atoms between adjacent lay- 
ers, but any matrix element which gives little weight to 
states along the nodal direction should give qualitatively 
similar results. We saw that the reduction in the den- 
sity of states caused by the pseudogap resulted in the 
low frequency region of the normal state optical conduc- 
tivity decreasing as temperature decreases instead of in- 
creasing as observed in-plane. This is in agreement with 
the experimental findings. This decrease continues in 
the superconducting state with the spectral weight redis- 
tributed to higher energies set by the pseudogap energy 
scale. This redistribution could serve as an explanation 
for the observation of the broad peak at 400cm^^ in un- 
derdoped YBCO (YBa2Cu306.7) which is not present in 
the optimally doped sample (YBa2Cu306.95).^^ 

We were also able to extract the behavior of the su- 
perfluid density from both a sum rule, and from the 
microwave conductivity. The superfiuid density behaves 
very differently for in-plane and out of plane, which man- 
ifested itself in the different shapes of the microwave con- 
ductivity as a function of temperature. The in-plane mi- 
crowave conductivity exhibits a peak, while the c-axis 
does not. A previous interpretation^^ took this obser- 
vation to mean that the c-axis interplane transport was 
incoherent. Here we attribute this to the very differ- 
ent temperature law associated with the c-axis super- 
fluid density as compared to the ab-plane. The shape 
of the c-axis conductivity constrains the temperature de- 
pendence of the scattering in the superconducting state. 
We also showed that the pseudogap suppresses the mi- 
crowave conductivity at high temperatures. It would be 
interesting to see more c-axis measurements, as they dis- 
play pseudogap physics more strongly than the in-plane 
counterparts. 
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